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Abstract 

A  mathematical  model  describing  the  distribution  of  the  concentrations  of  electrons  and  electron  holes  and  the  potential  along  the  thickness 
of  electrolyte  is  developed.  This  model  facilitates  the  study  of  various  electrolyte  materials  and  the  effects  of  operating  temperature,  oxygen 
partial  pressure  and  electrolyte  thickness  on  oxygen  semi-permeability.  The  model  is  applied  to  a  study  of  the  performance  of  a  ceramic 
oxygen  generators  (CoGs)  and  a  solid  oxide  fuel  cells  (SoFCs).  In  a  CoG,  CO  +  C02lelectroloytelpure  02  and  N2  +  02leleetroloytelpure  02 
are  considered  for  a  system  which  requires  low  02  and  high  02  semi-permeability,  respectively.  In  a  SoFC,  for  doped-Zr02  at  800  °C,  the 
thinner  the  electrolyte,  the  higher  the  energy  efficiency  and  the  output  power  density  will  be  at  intermediate  current  density  due  to  the  higher 
Ohmic  loss  of  the  electrolyte.  Thus,  obtaining  low  02  semi-permeability  for  a  thin  electrolyte  is  also  desirable  in  SoFC  development. 
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1.  Introduction 

Doped-Zr02,  doped-Ce02  and  doped-Bi203  are  ionic 
conductors  that  have  been  widely  considered  as  ideal  elec¬ 
trolyte  materials  for  ceramic  oxygen  generators  and  solid 
oxide  fuel  cells  (SoFCs)  [1-4].  Such  application  mainly 
considers  the  ionic  behaviour  of  the  electrolyte  materials. 
Though  they  are  ionic  solids,  the  diffusion  of  electronic 
components  (electrons  and  electron  holes)  across  the  elec¬ 
trolyte  cannot  be  neglected.  The  diffusivity  of  the  charge 
carriers  depends  largely  on  the  mobility  and  the  operating 
temperature.  Generally,  the  electronic  conductivity  is  con¬ 
sidered  to  arise  from  the  following  reactions  [5]. 


At  low  oxygen  pressures: 

>o  =  +2  +  V6  +  2e- 

0) 

At  medium  and  high  oxygen  pressures: 

02  +  Vq  =  Oq  +  2h 

(2) 

In  a  H2-02  fuel  cell  environment: 

I2  +  Oq  =  Vq  +  H20  +  2e_ 

(3) 

Under  normal  CoG  and  SoFC  operating  conditions,  the 
concentrations  of  oxygen  ions  and  oxygen  vacancies  remain 
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almost  constant  [6].  Thus,  based  upon  Eqs.  (1)  and  (2),  the 
concentrations  of  free  electrons  and  electron  holes  at  the 
gaslelectrolyte  interface  can  be  written  as: 

ce-  =  ke  pQf  (4) 

ch  =  khPo2  (5) 

where  the  reaction  rate  constants  kt  and  are  temperature- 
dependent  parameters. 

Since  there  is  mixed  ionic  and  electronic  conducting  beha¬ 
viour  in  these  oxide  electrolytes,  the  electronic  current  can 
flow  through  the  electrolyte  even  at  the  open-circuit  condition. 
Thus,  the  terminal  voltage  is  somewhat  lower  than  the  the¬ 
oretical  voltage  of  the  cell  reaction  and  the  performance  of  the 
cell  decreases  with  the  increase  of  electronic  conductivity  [7]. 

Though  electronic  conductivity  of  solid  oxides  is  undesir¬ 
able  for  the  electrolyte  used  in  SoFCs,  such  mixed-conduc¬ 
tive  behaviour  has  been  extensively  researched  for  gas- 
separation  membranes,  which  can  work  without  the  need 
of  electrodes  and  externally  applied  current  [8-10].  In  this 
case,  applying  a  differential  oxygen  pressure  across  the 
electrolyte  can  establish  concentration  gradients  of  free 
electron  and  electron  holes  across  the  electrolyte.  Thus,  free 
electrons  or  electron  holes  diffuse  through  the  ceramic 
membrane  naturally,  and  the  oxide  ions,  which  migrate  in 
the  counter  direction,  will  compensate  the  electrical  charge. 
The  net  effect  is  a  flux  of  oxygen  through  the  membrane, 
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Nomenclature 

ce-  concentration  of  free  electron  (mol  cm-3) 
c ^  concentration  of  electron-hole  (mol  cm-3) 

Ce  dimensionless  concentration  of  free  electron 

dimensionless  concentration  of  electron-hole 
Df  diffusion  rate  of  species  i  (cm2  s-1) 

Eq-  activation  energy  of  free  electron  (eV) 

Eh  activation  energy  of  electron-hole  (eV) 

F  Faradic  constant,  96  486  C  mol-1 

AHyQ  activation  energy  of  oxygen  vacancy  (eV) 

I  current  density  (A  cm-2) 

J  mass  flux  (mol  cm-2  s-1) 

kQ-  free-electron  reaction  rate  constant 

kh  electron  hole  reaction  rate  constant 

k  Boltzman  constant,  1.38  x  10-23  J  K_1 

L  electrolyte  thickness  (cm) 

p  oxygen  partial  pressure  (atm) 

Ra  area  polarisation  resistance  of  anode  (Q  cm2) 

Rc  area  polarisation  resistance  of  cathode  (Q  cm2) 

Re  area  electronic  resistance  (Q  cm2) 

Rj  area  ionic  resistance  (Q  cm2)  or  homogeneous 

reaction  rate  (mol  cm-3  s-1) 

T  temperature  (K) 

ut  mobility  of  species  i 

V  gas  volume  (cm3) 

z  electron  number 

Greek  symbols 

4>  inner  potential  of  the  electrolyte  (V) 

o  electrical  conductivity  (S  cm-1) 

v  gas  flow  velocity  (cm  s_1) 

Subscripts 
A  anode 

C  cathode 

e“  free  electron 

h  electron-hole 

Vq  oxygen  vacancy 


which  is  the  so-called  ‘oxygen  electrochemical  semi-perme¬ 
ability’  . 

In  this  study,  an  electrolyte  conduction  model  which 
considers  the  electronic  conductivity  is  developed  to  inves¬ 
tigate  the  transport  phenomena  of  free  electrons  and  electron 
holes  through  the  electrolyte,  and  to  evaluate  systemically 
the  effect  of  electronic  conduction  on  the  performance  of 
CoGs  and  SoFCs. 

2.  Electrolyte  model  considering  the  electronic 
conductivity 

With  reference  to  the  schematic  diagram  shown  in  Fig.  1, 
a  one-dimensional  electrolyte  model  was  developed  based 
on  the  following  assumptions: 


f) 


Fig.  1.  Schematic  diagram  of  electrolyte  reactions  with  transport  of  charge 
species. 

(i)  the  concentration  of  oxygen  vacancies  in  the  electrolyte 
is  constant  and  is  independent  of  the  oxygen  partial 
pressure  in  the  range  of  10“ 18  to  10  atm; 

(ii)  the  temperature  is  uniform  throughout  the  electrolyte; 

(iii)  the  ionic  and  electronic  conductivities  are  uniform 
throughout  the  electrolyte,  and  are  temperature- 
dependent  parameters; 

(iv)  the  current  and  species  concentrations  (electrons  and 
electron  holes)  vary  only  along  the  y-axis; 

(v)  the  mobility  of  the  charge  species  is  uniform 
throughout  the  electrolyte; 

(vi)  oxygen  surface  exchange  imposes  no  constraints  on 
the  ionic  transport. 

The  Nernst-Planck  equation  is  used  to  describe  the  flux  of 
species,  i.e. 

Ji  =  -ZiUiFci\/(j)  -  DiVci  (6) 

which  states  that  a  dissolved  species  can  move  by  migration 
and  diffusion.  The  current  density,  I  arising  from  the  motion 
of  the  charged  particles  is  given  by: 

/  =  (7) 

i 

All  charge  species  within  the  electrolyte  must  obey  the 
conservation  of  mass,  i.e. 

VJi+Ri  =  0  (8) 

Combining  Eqs.  (4)  and  (5)  can  eliminate  the  dependency  on 
oxygen  pressure,  thus 

ce-  X  rh  =  ks-pf2/4  x  khp'f  =  ke-kh  =  K(T )  (9) 

The  dimensionless  concentrations  of  free  electrons  and 
electron  holes  are  defined  as: 

Ce_f|y  and  Ch3f  (10) 

Ce-  C ^ 

The  requirement  of  electroneutrality  of  the  electrolyte 

yields: 

yy,  <•,  =  o  (ii) 

i 
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Since  the  concentration  of  oxygen  vacancies  is  constant 
(Vcy6  =  0and7?v5  =  0),  combining  Eqs.  (6)  and  (8)  gives 
Eq.  (12).  Assuming  homogenous  reaction  between  the  free 
electrons  and  electron  holes  within  the  electrolyte 
(Rq-  =  R^)  and  using  Nemst-Einstein  relation  Dt  =  UiRuT 
and  combining  Eqs.  (6),  (8)  and  (10)  written  for  respective 
free  electrons  and  electron  holes  yields  Eq.  (13).  Finally, 
combining  Eqs.  (9)  and  (10)  yields  Eq.  (14): 


V2</>  =  0 


(12) 


AO  (2)  d  T/t^) 


RuTc 


(2) 


z2F2 


-V2Ce-  =0 


(13) 


(14) 


To  solve  for  <f> ,  Ce-  and  Ch  and  the  following  six  boundary 
conditions  are  required: 


^  =  0 
/  (2)\  1/4 

ci-  =  1 


7(t)  =  _/( 2)  =  _(o-(2)Ce 


'“Ch 


RUT 


00 


■  cv( 
(0 


J(2) 

— veil'  +  ^-vc; 


G} 


.(0 


Ze- 


Zt 


(15a) 

(15b) 

(15c) 


(15d) 


ef  = 

1 

(15e) 

/  (2)\  */4 

e?  - 

fe 

(15f) 

\Po2J 

Eq.  (15d)  was  obtained  by  combining  Eqs.  (6),  (7)  and  (10) 
written  for  respective  charge  species.  The  conductivity, 
Gi  =  zfF2UiCi  at  the  electrolyte  boundary  can  be  expressed 
as  a  function  of  temperature  and  oxygen  partial  pressure 
[11],  i.e. 

<7v6  =  <7°6exp(^)  (16) 


Cq- 


O  -1/4 

°e~Po2  eXP 


O  1/4 

=  %Po2  exP 


(17) 

(18) 


Table  1  shows  typical  electrical  conductivities  and  activation 
energies  used  in  this  study. 

The  differential  Eqs.  (12)— (14)  do  not  yield  any  analytical 
solution  but  can  be  solved  by  a  numerical  method.  Once  the 
distribution  of  species  concentration  is  determined,  the  ionic 
current  density  and  electronic  current  density  (both  con¬ 
tributions  of  free  electrons  and  electron  holes)  flowing 
through  the  electrolyte  can  be  obtained.  Thus, 


Ii  = 

1 e  =  (f-}Ce-  +  AA)V4>(1) 


42) 


Ai) 


—  VC_;+^-VC 


A) 


Zq~ 


(19) 


(20) 


Under  open-circuit  conditions  ionic  current  is  equal  to  the 
electronic  current. 


3.  Results  and  discussion 

All  simulation  results  refer  to  the  operating  conditions  of 
SoFCs  andCoGs.  For  SoFCs,  the  oxygen  partial  pressure  is  set 
to  0.21  atm  at  the  cathode  side  and  10  18  atm  (balanced  by 
hydrogen  and  water  vapour  pressures)  at  the  anode  side.  For 
CoG,  two  systems,  CO  +  C02lelectrolytelpure  02  and  N2+ 
02lelectrolytelpure  02  are  considered.  The  former  is  an  indi¬ 
cative  gas  composition  for  Mars,  and  the  latter  for  Earth.  Two 
electrolytes,  doped-Zr02  and  doped-Bi203,  are  examined. 

The  distributions  of  the  dimensionless  concentration  of  free 
electrons  and  electron  holes  across  the  electrolyte  are  given  in 
Fig.  2.  On  the  high  oxygen  pressure  side  (y/L  =  0),  the 
conductivity  is  due  mainly  to  the  electron  holes,  while  on 
the  low  oxygen  pressure  side  (y/L  =  1 )  it  is  due  mainly  to  the 
free  electrons.  Since  an  equilibrium  condition  exists  between 
free  electrons  in  the  conduction  band  and  electron  holes  in  the 
valence  band  (annihilation),  the  intersection  of  the  concentra¬ 
tion  curves  of  free  electrons  and  electron  holes  is  expected  at  a 
certain  point,  3,  in  the  electrolyte.  This  critical  point  shifts 
from  the  high  (y/L  =  0)  to  the  low  oxygen  partial  pressure 
side  (y/L  =  1)  whenever or increases.  With  a  known 
distribution  of  electronic  charge  concentrations,  the  electro¬ 
nic  currents  due  to  the  respective  diffusion  of  free  electrons 


Table  1 

Electrical  conductivities  and  activation  energies  of  selected  electrolytes 


Electrolyte 

Temperature  range  (°C) 

0y  (S  cm  l) 

A Hy0  (eV) 

(Tg_  (S  cm  x) 

Ee-  (eV) 

(S  cm  x) 

(eV) 

Reference 

Zr02  +  8  mol%  Y203 

700-1200 

1.63  x  102 

0.79 

1.31  x  107 

3.88 

2.5 

1.67 

[12] 

BiO3+20  mol%  M203a 

500-700 

8.68  x  104 

1.08 

2.59  x  105 

2.24 

580.4 

1.12 

[13] 

M2O3:  Er203  or  Y2O3. 
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Fig.  2.  Distributions  of  dimensionless  concentration  of  free  electrons  and 
electron  holes  along  the  electrolyte  thickness. 

and  electron  holes  can  be  calculated;  the  results  are  given  in 
Fig.  3.  Along  the  thickness  of  electrolyte  from  the  high  to  the 
low  oxygen  partial  pressure,  the  electrolyte  changes  from  a  p- 
type  semiconductor  to  an  n-type  semiconductor.  This  phe¬ 
nomenon  can  be  used  to  distinguish  the  contribution  of  free 
electrons  and  electron  holes  to  the  total  electronic  conductiv¬ 
ity  of  the  electrolyte  by  semi-permeability  measurements  [5] . 

As  was  discussed  above,  the  electronic  current  due  to  the 
diffusion  of  free  electrons  and  electron  holes  is  compensated 
by  the  counter-flow  oxygen  ions.  The  rate  of  oxygen  semi¬ 
permeation  is  determined  by  the  electronic  current  that  flows 
through  the  electrolyte  and  can  be  expressed  as: 


where  /e  can  be  determined  from  Eq.  (20). 


The  effect  of  electrolyte  thickness  on  the  oxygen  semi¬ 
permeability  for  doped-Zr02  and  doped-Bi203  electrolytes 
is  presented  in  Fig.  4.  It  was  found  that  for  both  the  doped- 
Zr02  and  doped-Bi203,  the  oxygen  semi-permeability  is 
linearly  proportional  to  the  reciprocal  of  the  electrolyte 
thickness.  When  L  — >  oo,  the  oxygen  semi-permeability 
trends  to  zero.  Compared  with  the  doped-Zr02  at 
1000  °C,  the  oxygen  semi-permeability  of  doped-Bi203  at 
600  °C  for  a  0.2  mm  thickness  (1/L  =  5)  is  about  twice  as 
high.  It  is  also  noted  that  the  higher  the  operating  tempera¬ 
ture,  the  higher  is  the  semi-permeability.  As  we  will  discuss 
later,  the  semi-permeability  will  cause  different  performance 
characteristics  of  the  SoFC  and  CoG. 

To  study  the  impact  of  conductivity  on  oxygen  semi¬ 
permeability,  the  electronic  and  ionic  conductivities  of  the 
electrolyte  are  varied  arbitrary.  The  effect  of  the  conductivity 
on  the  oxygen  semi-permeability  using  the  values  from 
Table  2  is  given  in  Fig.  5. 

Comparing  curves  A  and  B,  though  the  ionic  conduc¬ 
tivity  is  increased  by  10  times  in  curve  B,  the  oxygen  semi¬ 
permeability  is  only  marginally  higher.  Comparing  curves 
A  and  C,  however,  if  both  the  free-electron  and  electron- 
hole  conductivities  are  increased  by  two  times  in  the  latter, 
the  corresponding  oxygen  semi-permeability  also  inc¬ 
reased  by  as  much.  This  phenomenon  is  particularly  useful 
for  the  selection  of  electrolyte  materials  for  use  in  gas- 
separation  devices.  For  ionic  domain  materials,  the  effec¬ 
tive  way  to  increase  the  oxygen  semi-permeability  is  to 
improve  the  electronic  conductivity  of  the  materials.  Some 
researchers  have  successfully  developed  high  oxygen  per¬ 
meation  composites  with  mixed  electronic  and  ionic  con¬ 
ductors  [14-16]. 

We  arbitrarily  assume  that  pure  C02  at  1  atm  is  fed  to  the 
cathode  side  of  a  CoG  simulating  the  Mars  atmosphere.  If 


Fig.  3.  Distributions  of  electronic  currents  along  electrolyte  thickness. 
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Table  2 

Electrical  conductivities  used  for  in  Fig.  5 


Curve 

cre-  (S  cm  x) 

<*h  (S  Cm  ^ 

(7y6  (S  cm  !) 

Aa 

1.77  x  1CT4 

1.82  x  10“4 

0.12 

B 

1.77  x  1(T4 

1.82  x  10^4 

1.2 

C 

3.54  x  1CT4 

3.64  x  10“4 

0.12 

aThe  conductivities  used  in  curve  A  are  based  on  doped-Zr02  at 
1000  °C  and  is  considered  as  the  baseline  for  comparison. 


the  concentration  gradients  of  C02  and  CO  at  the  cathode 
surface  are  neglected,  the  pressure  of  CO  on  the  cathode 
surface  with  a  drawn  current  density  of  I  and  a  C02  flow 
velocity  of  v  is  given  by 


P  co 


ncoRuT  _  ncoJcoKT  _  IRUT 

v  -  v  “TiT 


(22) 


Thus,  the  oxygen  partial  pressure  at  the  cathode  side  can 
be  calculated  for  any  temperature  based  on  equilibrium 


thermodynamics  [3],  i.e. 


P  o2 


x  1.75  x  104  x  exp 

Pco 


2.82  x  105\12 

rJt  J 


(23) 


At  open  circuit  (I  =  0)  po2  =  lpco2  and  Eq.  (23)  can  be 
rewritten  as: 


Po2  = 


Pc  o2 


x  3.5  x  104  x  exp  - 


2.82  x  105 


R„T 


2/3 


(24) 


The  effect  of  temperature  on  the  oxygen  semi-permeability 

for  different  electrolyte  thicknesses  at  open-circuit  is  shown 

in  Fig.  6a  and  b.  In  Fig.  6a,  is  calculated  from  Eq.  (24) 

and  p is  set  to  1  atm.  To  simulate  the  CO  +  C02lelectro- 
u2  (2) 

CoG  system,  while  in  Fig.  6b py0 ^  and  jPq2  are  set  to  0.21  and 

1  atm,  respectively,  to  simulate  the  N2  +  02 1  electrolyte  I  pure 
02  CoG  system.  Results  show  that  the  rate  of  oxygen 
permeation  of  doped-Bi203  electrolyte  with  a  thickness  of 
0.2  mm  increases  monotonically  with  increasing  tempera¬ 
ture  and  reaches  1.523  x  10-8  and  0.202  x  10-8  mol 


Fig.  5.  Sensitivity  of  electrical  conductivity  on  oxygen  semi-permeability. 


S.H.  Chan  et  al.  /  Journal  of  Power  Sources  111  (2002)  320-328 


325 


450  500  550  600  650  700  750 


(a) 


Temperature  (°C) 


(b)  Temperature  (°C) 

Fig.  6.  (a)  Effect  of  temperature  on  oxygen  permeation  at  open-circuit  for  CO  +  C02ldoped-Bi203l02  CoG  system,  (b)  Effect  of  temperature  on  oxygen 
permeation  at  open-circuit  for  N2  +  02ldoped-Bi203l02  CoG  system. 


cm2  s-\  respectively,  at  700  °C  for  the  two  systems.  For  all 
temperatures,  more  oxygen  permeation  can  be  expected  for 
thinner  electrolytes,  which  agrees  with  the  results  shown  in 
Fig.  4. 

The  effect  of  current  density  on  the  rate  of  oxygen 
generation  for  the  CO  +  C02lelectrolytelpure  02  CoG  sys¬ 
tem  is  given  in  Fig.  7.  Due  to  the  oxygen  semi-permeability 
from  the  pure  oxygen  side  to  the  CO  +  C02  side,  which  is  in 
the  opposite  direction  of  the  ionic  flow  due  to  current 
applied,  the  actual  oxygen  produced  in  doped-Bi203  at 
700  °C  is  lower  than  the  amount  expected.  Because  of  the 
low  electronic  conductivity  in  doped-Zr02  at  1000  °C,  the 
actual  oxygen  produced  is  almost  same  as  the  theoretical 
value.  In  the  case  of  the  N2  +  02 1  electrolyte  Ipure  02  CoG 
system,  the  migration  of  oxygen  ions  due  to  semi-permea¬ 
tion  is  in  the  same  direction  as  the  migration  of  oxygen  ions 
due  to  current  applied.  Hence,  the  oxygen  semi-permeability 
behaviour  of  the  electrolyte  is  actually  beneficial  to  the 
generation  of  the  oxygen.  Results  in  Fig.  8  show  that 
the  doped-Bi203  based  CoG  can  actually  produce  more 


oxygen  than  the  theoretical  amount.  Thus,  selection  of 
electrolyte  materials  for  CoG  application  should  consider 
the  operating  environment  and  the  working  principle  of  the 
CoG  to  maximise  the  oxygen  yield. 

When  ion-conducting  membrane  is  used  in  a  SoFC,  efficient 
operation  of  the  fuel  cell  requires  the  ionic  conductivity  to  be 
much  higher  than  its  electronic  counterpart  in  the  electrolyte. 
The  equivalent  circuit  of  a  SoFC  is  given  in  Fig.  9  and  includes 
consideration  of  electronic  conduction  in  the  electrolyte.  The 
output  power  density  of  the  fuel  cell  is  given  by 

Pout  —  /out  X  £out  (25) 

The  rate  of  total  energy  input  to  the  SoFC  is  the  summa¬ 
tion  of  the  electrical  power  output  (PG ut)  and  the  Joule 
heating  due  to  the  Ohmic  losses  associated  with  ionic  flow 
(Pi),  electronic  flow  (Pe),  anode  polarisation  (PA)  and  cath¬ 
ode  polarisation  (Pc).  Thus, 

P  total  =  Pi+Pe+^A+Pc+^>out 

=  fout^out  +  lfR{  +  /2Pe  +  l\R  a  +  I2cRc 


(26) 
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Fig.  7.  Effect  of  current  density  on  oxygen  generation  rate  for  different  electrolyte  materials  used  in  CO  +  C02ldoped-Bi203l02  CoG  system. 


Fig.  8.  Effect  of  current  density  on  oxygen  generation  rate  for  different  electrolyte  materials  used  in  N2  +  C02ldoped-Bi203l02  CoG  system. 


Anode  Cell  Cathode 


Fig.  9.  Equivalent  circuit  of  SoFC  including  electronic  conduction  of 
electrolyte. 


The  energy  efficiency,  s ,  is  defined  the  ratio  of  the  electrical 
power  output  to  the  rate  of  total  energy  input  to  the  fuel  cell, 
i.e. 


j nut  Erf 


8  = 


total 


IautEa 


/^1+/e27?e+/^A+/^c 


(27) 


The  effect  of  the  current  density  on  energy  efficiency  and 
output  power  density  is  illustrated  in  Figs.  10  and  11, 
respectively,  for  doped-Zr02  at  800  °C.  The  thinner  the 
electrolyte,  the  higher  the  output  power  density  and  energy 
efficiency  will  be  at  intermediate  current  density.  In  addi¬ 
tion,  the  peak  energy  efficiency  shifts  towards  lower  current 
densities  for  a  thicker  electrolyte.  As  shown  above  (Fig.  4), 
the  rate  of  oxygen  semi-permeability  (or  the  electro¬ 
nic  current  density)  decreases  with  increase  in  electrolyte 
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Current  Density  (A/cm2) 

Fig.  10.  Energy  efficiency  as  a  function  of  current  density  for  different  doped-Zr02  electrolyte  thickness  (operation  at  800  °C). 


Fig.  11.  Output  power  density  as  a  function  of  current  density  for  different  doped-Zr02  electrolyte  thickness  (operation  at  800  °C). 


thickness.  Thus,  the  effect  of  electrolyte  thickness  on  energy 
efficiency  depends  on  the  combined  losses  due  to  the 
electrolyte  resistance  and  the  oxygen  semi-permeability. 
For  an  electrolyte  of  ionic  conductivity,  such  as  doped- 
Zr02,  the  Ohmic  loss  dominates  the  energy  efficiency  of 
the  fuel  cell  at  intermediate  current  density. 

4.  Conclusions 

An  electrolyte  model  for  ion-conducting  membranes  used 
in  ceramic  oxygen  generators  and  solid  oxide  fuel  cells  has 
been  developed.  The  model,  which  consists  of  a  set  of  three 
second-order  differential  equations  and  six  boundary  con¬ 
ditions,  establishes  the  relationships  between  the  potential 
and  the  concentrations  of  free  electrons  and  electron  holes 
along  the  thickness  of  the  electrolyte.  It  facilitates  the  study 
of  various  electrolyte  materials  as  well  as  the  effects  of 


operating  temperature,  oxygen  partial  pressure  and  electro¬ 
lyte  thickness  on  the  oxygen  semi-permeability.  The  model 
is  used  to  analyse  the  effect  of  oxygen  semi-permeability  on 
the  performance  of  a  CoG  and  SoFC.  The  contribution  of 
free  electrons  and  electron  holes  to  the  total  electronic 
conductivity  depends  highly  on  the  oxygen  partial  pressure 
at  the  gas  I  electrolyte  interface.  For  ionic  domain  electrolyte, 
the  oxygen  semi-permeability  is  mainly  influenced  by  the 
electronic  conductivity  of  the  electrolyte.  An  effective  way 
to  enhance  the  throughput  of  a  differential  oxygen  generator 
is  to  use  a  material  with  high  electronic  conductivity  for 
improved  oxygen  semi-permeability.  In  addition,  high  oxy¬ 
gen  semi-permeability  can  be  achieved  with  high  operating 
temperature,  high  differential  oxygen  pressure,  or  a  thinner 
electrolyte.  For  a  CO  +  C02lelectrolytelpure  02  CoG  sys¬ 
tem,  it  is  desirable  to  use  an  electrolyte  with  low  oxygen 
semi-permeability  because  the  oxygen  semi-permeation  is 
in  counter  direction  to  the  ionic  flow  under  applied  current. 


328 


S.H.  Chan  et  al./  Journal  of  Power  Sources  111  (2002)  320-328 


For  a  N2  +  02lelectrolytelpure  02  CoG  system,  however, 
high  oxygen  semi-permeability  is  desirable  as  oxygen  semi¬ 
permeation  is  in  the  same  direction  of  the  ionic  flow  under 
applied  current.  In  a  SoFC  with  doped-Zr02  at  800  °C,  the 
thinner  the  electrolyte,  the  higher  is  the  energy  efficiency  at 
intermediate  current  density  because  of  the  higher  Ohmic 
loss  of  the  electrolyte. 
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